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radiation hazard from backflow of fission fragments 


FROM THE PLUME OF A GAS-CORE NUCLEAR ROCKET 
by CharJes C. Masser 
Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio 

SUMMARY 

Analytical calculations are performed to determine the particle flux 
that leaves a rocket plume in a vacuum and strikes the rocket forward of 
the nozzle. The rocket plume is generated by the products of a gas-core 
nuclear reactor and consists of hydrogen, uranium, and fission fragments. 
The backflow flux is determined by the number of particles that leave the 
continuum of the plume and strike the rocket without experiencing any par- 
ticle collisions. A percentage of these particles that strike the rocket are 
fission fragments and the radiation dose to the crew is calculated. 

A total of 1. 68 pounds of fission fragments are formed from a gas- 
core rocket that produces one million pounds of thrust at a specific im- 
pulse of 1500 seconds for a propellant consumption of one million pounds. 

O 

Of the 1. 68 pounds of fission fragments only 0. 42x10’ pound flow back 

from the plume. And only 0.66x10”*® pound of fission fragments strike 

the rocket. This weight of fission fragments causes a total integrated 

-4 

dose of less than 1. 3x10 rad to the crew, which is well below acceptable 
levels. 


INTRODUCTION 

In the open cycle gas-core nuclear rocket concept (Fig. 1) the heat 
source is fissioning uranium gas. This released heat is radiated to and 
absorbed by the hydrogen propellant. The heated propellant is exhausted 
through a nozzle, producing thrust. In an open cycle gas-core rocket the 
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fission fragments that are formed and the unfissioned uranium fuel are 
.ilso exhausted into the vacuum of space. As the exhaust plume is formed, 
.i small percentage of the plume products have a sufficient velocity and the 
proper direction to leave the plume and flow back towards the rocket. It 
is the purpose of this paper to evaluate the radiation hazard to the crew 
associated with a backflow of fission fragments that strike the rocket. It 
is assumed that the fission fragments striking the rocket stick to it. 

There are other radiation sources associated with the gas-core reac- 
tor such as radiation from the reactor core and radiation from inside the 
plume itself. These radiation sources along with solar radiation must be 
ultimately considered. This study, however, is concerned only with that 
part of the total radiation problem that arises from fission fragments 
leaving the plume and striking the rocket. 


ROCKET ENGINE CHARACTERISTICS 

Certain exit nozzle conditions are selected in order to make a specific 
calculation of the backflow of fission fragments from the plume. Roback 
(Ref. 1) calculated performance parameters for hydrogen at various stag- 
nation pressures and temperatures. The selection was made to pick per- 
formance parameters which produced a high flow of fission fragments, 
thereby maximizing the radiation problems. A one million pound propel- 
lant storage capacity was picked, which is typical of a manned Mars mis- 
sion. 

Thrust, lb 

Specific impulse, sec . . 

Chamber pressure, atm . 

Chamber temperature, °R 
Reactor power, MW . . . 

From table 168 of Roback (Ref. 1) for a 

Pressure nozzle exit. _ ^q-3. 

Pressure nozzle throat 


1 000 000 
. . .1500 
. . .1000 
. 10 000 
. 41 900 
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Kxit temperature, °R % 3460 

Exit density, lb/ft 3 7.978xl0~ 4 

Kxit velocity, ft/sec 48 507 

Exit Mach number 4.587 

Exit molecular weight . 2.015 

Molecular diameter, cm 2.75xl0“ 3 


For the given thrust of one million pounds and a specific impulse of 
1500 seconds, the propellant flow rate is 667 pounds per second. There- 
fore, if the propellant storage for the mission is one million pounds, the 
total firing time is 1500 seconds. 


CALCULATION OF PLUME DENSITY 


It has been shown that in nozzle plume flows the mass flux pV 
varies inversely as the square of the distance from the source point. 
Hill and Draper (Ref. 2) have shown that the density in the plume can be 
closely approximated by 


p = 


4p MB 
e e 



e^ e -A^(l-cos 0Y 
(y+l)/2(y-l) \rj 


(1) 


where the coordinates of r and 0 are shown in Fig. 2; p is density 
at any poi: ; in the plume; p g is exit density; is exit Mach number; 
y is the ratio of specific heats; and B and X are constants. 

Also from Hill and Draper (Ref. 2) we have 

i/(y-D 



( 2 ) 
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Where Cp and 
Eqs. (4. 33) and 


C F are thrust coefficients and are evaluated using 
r max 

(4.34) of Shapiro (Ref. 3). 


DEFINITION OF PLUME BOUNDARY 

The gas inside the plume is considered to be a continuum. As the 
plume expands, the density decreases to a point where free molecular 
flow can be assumed. At this point we can define our plume boundary. 
From Grier (Ref. 4), the plume boundary is defined by a surface on 
which the Knudsen number is constant. The Knudsen number, K, is de- 
fined as 

K = — (4) 

r 

where r is the distance from the plume source point to the plume bound- 
ary, and l is the mean free path of the particle and is given by Santeler 
et al. (Ref . 5) as 


l = 3.736X10' 25 — (5) 

P6 2 

where m is molecular weight ar.d 5 is molecular diameter. 

Solving Eq. (5) for density and substituting in Eq. (4) we have 

p = 3.736X10" 25 (6) 

K6 2 r 

Combining Eq. (6) and Eq. (1) we have the equation of a constant Knudsen 
surface, which defines the boundary of the*plume. 
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r 

r„ 


1.07xl0 25 p M BK5 2 r 
e c e 


-X 2 (l-cos 0) 2 


m 1 + 


•211-1 M 2 


1/2 


(y+l)/2(y-l) 


(7) 


<r + h 


In this equation the maximum 6 was taken to be the maximum 
Prandtl-Meyer expansion angle for a Mach number of 4, 87. The backflow 
of flux is determined by the number of particles that leave this defined 
boundary and strike a particular point P(x, y, z). A model of this is 
shown in Figs. 2 and 3. 


CALCULATION OF BACKFLOW FLUX 


Noller (Ref. 6) has calculated the density at a point P(x, y, z) due to 
the molecules leaving an elemental volume dr adjacent to a constant 
Knudsen surface. Noller ’s result was: 

^ e- u2 /l U ,os * + fl + U 2 ccs 2 ^ 
n T \2 \2 

Xe U 2 cos 2 ^AP\ 


[1 + erf(U cos i//)]\d£2 


( 8 ) 


where n T is the number density of molecules in dr, \}/ is the angle be- 
tween the direction of V and the direction of flight of the molecules con- 
tributing to the density at P(x, y,z), U = v/(2RT/m), R is the gas con- 
stant, T is the temperature, and df2 is the differential solid angle of 
dr as seen from P(x, y, z). Figure 3 shows the coordinate system used. 

Equation (8) was recalculated in this study to express mass flux (F) 
instead of density. The following equation resulted: 




-2/3 2 / 

ZZ = e"^ \ (1 + U 2 cos 2 i//) + a[v U cos ip 

"T 2 l 

+ U 2 cos 2 t//j[l + erf(l + cos i//)] (9) 


From Grier (Ref. 4), d£2 is given by 


dft = 


-2- sin 6 cos cp( 1 + 2A^ cos 6 - 2A 2 cos 2 0) 


r 

r 


,2 2 „ ^2 . 2 , 


\2 /r \ 3 

+ — (cos 0 - 2A* sin*0 + 2A* sin*# cos 6} f— ] (-^] sin 0 d# dtp 
r„ ' r e/ \ L / 


( 10 ) 


Where L is the vector from dr to the point P(x,y,z). Also, from 
Grier (Ref. 4) |l| is given by 


— - = rliL sin 0 cos (p- — + — cos 6 1 


r e \ r e 


r r « 
e e 


+ 6\%- cos 0 cos <r - ~ sin 0 ) * <p ~ sin (p (11) 


The particle mass flux at any point upstream of the nozzle face is 

(p=(n/2) 


-3/2 tt 2 / t6M max 
F - - e‘ L 



n /() + U 2 cos 2 i//) + aJti U cos i// 


'0=0 w ^=-(r/2) 


+ U 2 cos 2 !/^ [1 + erf(U cos p) je u2 cos2 ^ 


dft (12) 


CALCULATION OF FISSION FRAGMENT FORMATION 


The total engine running time is calculated by dividing the mass flow 
of propellant by the total mass of propellant carried. This results in a 
total engine running time of 1500 seconds. The number of fission frag- 
ments formed is calculated using Eq. (2.49) of Glasstone and Sesonski 
(Ref 7) 


Reactor power (watts) - E e - r g e - c0 - n j (13) 

3.1X10 10 

Since reactor power and engine running time are known, the number 
of fission fragments are known. It is also assumed the average molecular 
weight of the fission fragments is 117.5. Therefore we have the values 
shown in the following table. 


■ 

Number of particles 

Weight cf particles, 
lb 

Fission fragments 
Hydrogen propellant 
Unfissioned uranium 

3.896X10 24 

1.36X10 32 

1.163X10 28 

1.68 

10 6 

10 4 


It is assumed the fission fragments are uniformly distributed through- 
out the plume which is composed chiefly of hydrogen, and that the number 
of fission fragments flowing from the plume is of the same relative con- 
centration as the plume This assumption should give an overestimate of 
the number of fission fragments leaving the plume because of the higher 
mobility of the hydrogen molecules. 


MASS FLUX TO ROCKET 

In order to calculate the number of fission fragments that strike and 
stick to the rocket, a particular shape was pic 1 eU. The configuration 
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picked was a cylinder 33 feet in diameter and 400 feet long. This size is 
sufficient for a tankage volume needed to store one million pounds of hy- 
drogen whicn is typical of a manned Mars mission. The crew is assumed 
to be located at the front of the rocket 400 feet from the nozzle. 

The distribution of fission fragments over the rocket is computed so 
the distance from the fission fragments to the crew is known. In Fig. 4 
we have an enlarged view of the rocket nozzle area. The plume boundary 
is for the case of Knudsen number equal to 1. For this case the total 

number of fission fragments formed is 1.68 pounds. Of this amount only 

-8 ft 

0.42x10 pound flow back across the plane z = 0. Of this 0.42x10 

pound, 0. 66x10" ^ pound deposits on the end surface of the rocket. It 

can easily be seen that only a small percentage of the total backflow ever 

strikes the rocket. 


RADIATION DOSE TO THE CREW 

The radiation dose to the crew is calculated in two parts: the radia- 
tion from the fission fragments that strike the end surface and the radia- 
tion from the fission fragments that strike the side of the rocket. In all 
cases studied the radiation dose level to the crew from the fission frag- 
ments striking the side of the rocket is at least six orders of magnitude 
less than the dose from the fission fragments striking the end of the 
rocket. Therefore the radiation from the fission fragments striking the 
end of the rocket are only discussed. 

In calculating the radiation dose several assumptions are used in 
order to determine the highest dose that cne could expect from the back- 
flow. First, the gamma radiation from the fission fragments is released 
after they fiow back and strike the rocket. Second, the fission fragments 
of the uranium atom release their total 6 MeV of energy (Glasstone and 
Sesonske, Ref. 7) in decay in 1 second. Third, an unshielded point ker- 
nel calculation was used to calculate the dose to the crew. 

The total radiation dose to the crew for these extremely conservative 
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assumptions in the case where the Knudsen number K - 1 is 1.3x10 
rad. For a Knudsen number of 1. 0 the mean free path of the particles is 
equal to the distance from nozzle. This is shown on Fig. 4. This level 
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of rad-aiion is small so no hazard to the crew is expected. Since the 

amount of fission fragments that strike the rocket may depend on the 

definition oi the plume boundary, a variation of this definition is shown 

on Fig. 5. The plume boundaries represent the cases where the Knudsen 

number K is equal to 0. 01, 0.1, and 1.0. The total radiation doses for 

-4 -4 -4 

these three cases are 8.3x10 , 3.6x10 , and 1.3x10 rad, resoec- 

tively. However, even for the case where K = 0.01 the radiation level 
*o the crew is small. Therefore, from the cases studied the radiation 
dose to the crew from fission fragments that flow back from the plume is 
considered not to be hazardous. 
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APPENDIX - SYMBOLS 

defined by Eq. (2) 
thrust coefficient 
maximum thrust coefficient 

mass flux 
Knudsen number 

vector distance from dr on constant Knudsen surface to 
point P(x,y,s) 

mean free path of gas molecules 
Mach number 
molecular weight 
particle density at point P(x, y, z) 
particle density in volume element 
gas constant 
spherical coordinates 
temperature 
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u 

V 

x, y, * 

v 

* 

6 

A 


P 

* 

Subscript: 

e 


dimensionless velocity | V ly^RT/m)*^ 

mean velocity of particles in plume 

Cartesian coordinates 

ratio of specific heats 

molecular diameter of molecules 

defined by Eq. (3) 

mass density 

angle between direction of V and direction 
nozzle exit 


of 


L 
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KNUDSFN TOTAL BACKFLOW OF FISSION FRAGMENTS RADIATION OOSC 
NUMBER FISSION FRAGMENTS ON ROCKET TO CfE.V 


0. 42x10"® LB 
1.16x10"® LB 
2.69x10 ® LB 


0. 66x10" L B 
1.83x]0' 10 LB 
4. 22xlO" 10 lB 


1 3*10^ PAD 
3 6»!0"‘ 5 RAO 
8 3 x10' 4 fi-'-D 


Figure 5. - Distribution of fission fragments for a plume Knudsen number of 1.0, 0. 1, am 0.01. 



